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Experimental and theoretical progress concerning the rare decay tt — > e + e~ is briefly reviewed. It includes the 
latest data from KTeV and a new model independent estimate of the decay branching which show the deviation 
between experiment and theory at the level of 3.3a. 



1. Introduction 

Astrophysics observables tell us that 95% of the 
matter in the Universe is not described in terms 
of the Standard Model (SM) matter. Thus, the 
search for the traces of New Physics is a funda- 
mental problem of particle physics. There are two 
strategies to look for the effects of New Physics: 
experiments at high energy and experiments at 
low energy. In high-energy experiments it is con- 
sidered that due to a huge amount of energy the 
heavy degrees of freedom presumably characteris- 
tic of the SM extension sector are possible to ex- 
cite. In low-energy experiments it is huge statis- 
tics that compensates the lack of energy by mea- 
suring the rare processes characteristic of such ex- 
tensions. At present, there is no any evidence 
for deviation of SM predictions from the results 
of high-energy experiments and we are waiting 
for the LHC epoch. On the other hand, in low- 
energy experiments there are rough edges indi- 
cating such deviations. The most famous exam- 
ple is the muon (g — 2). Below it will be shown 
that due to recent experimental and theoretical 
progress the rare process ir° — > e + e~ became a 
good SM test process and that at the moment 
there is a discrepancy between the SM prediction 
and experiment at the level of 3.3cr deviation. 
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2. KTeV data 

In 2007, the KTeV collaboration published the 
result [T] for the branching ratio of the pion decay 
into an electron-positron pair 
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) = (7.48 ±0.38) • 10~ 8 . (1) 



The result is based on observation of 794 candi- 
date 7T° — > e+e~ events using Kl — > 37r° as a 
source of tagged 7r°s. Due to a complicated chain 
of the process and a good technique for final state 
resolution used by KTeV this is a process with low 
background. 



3. Classical theory of ir° 



decay 



The rare decay 7r° —> e + e~ has been stud- 
ied theoretically over the years, starting with the 
first prediction of the rate by Drell [2] . Since no 
spinless current coupling of quarks to leptons ex- 
ists, the decay is described in the lowest order of 
QED as a one-loop process via the two-photon 
intermediate state, as shown in Fig. 1. A fac- 
tor of 2 (m e /m T ) 2 corresponding to the approxi- 
mate helicity conservation of the interaction and 
two orders of a suppress the decay with respect 
to the 7T° — > 77 decay, leading to an expected 
branching ratio of about 10~ 7 . In the Standard 
Model contributions from the weak interaction to 
this process are many orders of magnitude smaller 
and can be neglected. 

To the lowest order in QED the normalized 
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One can attempt to reconstruct the full ampli- 
tude by using a once-subtracted dispersion rela- 
tion 5 



A(q 2 )=A(q 2 = 0) + 
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(6) 
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If one assumes that Eq. Q is valid for any q 2 , 
then one arrives for q 2 > 4m 2 at [6I7I8| 



Figure 1. Triangle diagram for the ir° — > e + e~ 
process with a pion 7r° — ► 7*7* form factor in the 
vertex. 



branching ratio is given by 
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0.988. The amplitude „4 can be written as 
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where q 2 = to 2 ,j> 2 = to 2 . i^ 7 * 7 » is the form 
factor of the transition 7r° — > 7*7* with off-shell 
photons. 

The imaginary part of A is defined uniquely as 

lmA ^ = w^F) ln{ve{q2))l (4) 



Pe{q 2 Y 

It comes from the contribution of real photons in 
the intermediate state and is model independent 
since F^*^» (0,0) = 1. Using inequality \A\ 2 > 
(Im.4) one can get the well-known unitary bound 
for the branching ratio 3 

B (ir° -> e + e~) (5) 

> ^unitary ^0 
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where Li 2 (z) = — f Q (dt/t) In (1 — t) is the dilog- 
arithm function. The second term in Eq. ([7]) takes 
into account a strong q 2 dependence of the am- 
plitude around the point q 2 — occurring due to 
the branch cut coming from the two-photon inter- 
mediate state. In the leading order in [me/m^) 2 , 
Eq. ([7]) reduces to 



Re.4 {ml) =A(q 2 = 0) 
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Thus, the amplitude is fully reconstructed up 
to a subtraction constant. Usually, this constant 
containing the nontrivial dynamics of the process 
is calculated within different models describing 
the form factor F n (k 2 ,q 2 ) [4151719110] . However, 
it has recently been shown in [10] that this con- 
stant may be expressed in terms of the inverse 
moment of the pion transition form factor given 
in symmetric kinematics of spacelike photons 



A (q 2 = 0) = 3 In 



(9) 
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Here, /i is an arbitrary (factorization) scale. One 
has to note that the logarithmic dependence of 
the first term on /i is compensated by the scale 
dependence of the integrals in the brackets. In 
this way two independent processes becomes re- 
lated. 



3 



4. Importance of CLEO data on -F,r 7 * 7 

In order to estimate the integral in Eq. ©, 
one needs to define the pion transition form fac- 
tor in symmetric kinematics for spacelike pho- 
ton momenta. Since it is unknown from the first 
principles, we will adapt the available experimen- 
tal data to perform such estimates. Let us first 
use the fact that F^y*^ (t,t) < F 7rT » 7 « (i,0) for 
t > in order to obtain the lower bound of 
the integral in Eq. ([9]). For this purpose, we 
take the experimental results from the CELLO 
[TT] and CLEO [T^] Collaborations for the pion 
transition form factor in asymmetric kinemat- 
ics for spacelike photon momentum which is well 
parametrized by the monopole form [12] 
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-23.2 ±0.1. 



Thus, for the branching ratio we are able to estab- 
lish the important lower bound which consider- 
ably improves the unitary bound given by Eq. ([5]) 
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It is natural to assume that the monopole form 
is also a good parametrization for the form factor 
in symmetric kinematics 
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The scale s% can be fixed from the relation for 
the slopes of the form factors in symmetric and 
asymmetric kinematics at low t [13J, 
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that gives s\ = so/2. Note that a similar reduc- 
tion of the scale is also predicted by OPE QCD 
from the large momentum behavior of the form 
factors: s? PE = s° PE /3 Q3]. Thus, the estimate 
for A (0) can be obtained from Eq. (jlip by shift- 
ing the lower bound by a positive number which 
belongs to the interval [31n(2)/2,31n(3)/2] 



A (<z 2 = 0) 



In 



si 



-21.9±0.3.(15) 



With this result the branching ratio becomes 
B (tt -> e+e-) = (6.23 ± 0.09) • 10" 8 . (16) 



This is 3.3 standard deviations lower than the 
KTeV result given by Eq. Q. 



Figure 2. Evolution of model predictions and 
comparison with the latest KTeV result. 



For this type of the form factor one finds from 
Eq. © that 



A(q 2 = 0) 



(11) 



5. Other decay modes 

The r\ — > decay can be analyzed in a sim- 
ilar manner. As in the pion case, the CLEO Col- 
laboration has parametrized the data for the r\- 
meson in the monopole form [T2]: 
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which is very close to the relevant pion parameter. 
Then following the previous case (with evident 
substitutions) , one finds the bounds for the q 2 — > 
limit of the amplitude i] — > (J- + fi~ as 

A v {q 2 = 0) (18) 
o / „CLEO \ r 

and for r] — » e + e~ one gets again Eq. (TTT]) . The 
obtained estimates allow one to find the bounds 
for the branching ratios 

B(r]^n + (i-) < (6.23 ±0.12)- 1CT 6 , (19) 
B (r? -> e+e - ) > (4.33 ± 0.02) • 10" 9 . 

It is important to note that for the decay r\ —> 
H + pT wc get the upper limit for the branching. 
This is because the real part of the amplitude for 
this process taken at the physical point q 2 = m 2 
for the parameter Sq^ eo remains negative and 
a positive shift due to the change of the scale 
sq v — > sir! reduces the absolute value of the real 
part of the amplitude |Re^4(m^)J. At the same 
time, considering the decays of n and r] into an 
electron-positron pair, the evolution to physical 
point ([7]) makes the real part of the amplitude 
to be positive for the parameter Sq LEO and the 
absolute value of the real part of the amplitude 
increases in changing the scales of the meson form 
factors. Thus, it would be very interesting to 
check experimentally the predicted bounds for the 
process r\ — ► 

The predictions for the decays rj — > l + l~ ob- 
tained by reducing the scale Sq^ eo — > si v for 
the case of the ry-meson transition form factor are 
given in Table [1] Also note the recent measure- 
ment by the WASA/Celcius Collaboration [15] 
which improves the upper limit for the branch- 
ing r\ — > e+e - . 

6. Possible explanations of the effect 

Therefore, it is extremely important to trace 
possible sources of the discrepancy between the 
KTeV experiment and theory. There are a few 
possibilities: (1) problems with (statistic) exper- 
iment procession, (2) inclusion of QED radiation 
corrections by KTeV is wrong, (3) unaccounted 



mass corrections are important, and (4) effects of 
new physics. At the moment, the last possibilities 
were reinvestigated. In [17] . the contribution of 
QED radiative corrections to the ir° — » e + e~ de- 
cay, which must be taken into account when com- 
paring the theoretical prediction (|16[) with the 
experimental result ([1]), was revised. Comparing 
with earlier calculations [18] , the main progress is 
in the detailed consideration of the 7*7* — > e + e~ 
subprocess and revealing of dynamics of large and 
small distances. Occasionally, this number agrees 
well with the earlier prediction based on calcula- 
tions [18] and, thus, the KTeV analysis of radia- 
tive corrections is confirmed. In [19 it was shown 
that the mass corrections are under control and 
do not resolve the problem. So our main conclu- 
sion is that the inclusion of radiative and mass 
corrections is unable to reduce the discrepancy 
between the theoretical prediction for the decay 
rate (fTl)]) and experimental result |T]). 

7. 7T° — > e + e decay as a filtering process 
for low mass dark matter 

If one thinks about an extension of the Stan- 
dard Model in terms of heavy, of an order of 100 
GeV or higher, particles, then the contribution of 
this sort of particles to the pion decay is negligi- 
ble. However, there is a class of models for de- 
scription of Dark Matter with a mass spectrum of 
particles of an order of 10 MeV [50]. This model 
postulates a neutral scalar dark matter particle 
X which annihilates to produce electron/positron 
pairs: \X ~ > e + e~. The excess positrons pro- 
duced in this annihilation reaction could be re- 
sponsible for the bright 511 keV line emanating 
from the center of the galaxy [3T]. The effects 
of low mass vector boson U appearing in such 
model of dark matter (Fig. [3|) were considered in 
22J where the excess of KTeV data over theory 
put the constraint on coupling which is consistent 
with that coming from the muon anomalous mag- 
netic moment and relic radiation [23]. Thus, the 
pion decay might be a filtering process for light 
dark matter particles. 

Further independent experiments at KLOE, 
NA48, WASAatCOSY, BES III and other facili- 
ties will be crucial for resolution of the problem. 
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Table 1 

Values of the branchings B (P — > l + l~) obtained in our approach and compared with the available ex- 
perimental results. 



B 


Unitary bound 


CLEO bound 


CLEO+OPE 


Experiment 


B e+e-) x 10 8 


> 4.69 


> 5.85 ±0.03 


6.23 ±0.09 


7.49 ±0.38 PQ 


B(r) -> M + M~) x lO 5 


> 4.36 


< 6.23 ±0.12 


5.11 ±0.20 


5.8 ±0.8 16] 


5(7? -» e+e - ) x 10 y 


> 1.78 


> 4.33 ±0.02 


4.60 ±0.06 


< 2.7 x 10 4 [15, 




Figure 3. Loop diagram for 7r° — > e + e process 
induced by the low mass exotic U* boson. 



Also important is to get more precise data on the 
pion transition form factor in asymmetric as well 
in symmetric kinematics. 
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